Proposal for an optical laser producing light at half the Josephson frequency 
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We describe a superconducting device capable of producing laser light in the visible range at 
half of the Josephson generation frequency with the optical phase of the light locked to the su- 
perconducting phase difference. It consists of two single-level quantum dots embedded in a p-n 
semiconducting heterostructure and surrounded by a cavity supporting a resonant optical mode. 
We study decoherence and spontaneous switching in the device. 

PACS numbers: 42.55.Px 73.40.-c, 74.45. +c, 78.67.-n 



Lasers and superconductors are both systems with 
macroscopic quantum coherence. In lasers, photons form 
a coherent state induced by stimulated emission of a 
driven system into a cavity mode. The resulting visi- 
ble coherent light is characterized by an optical phase 
1]. In superconductors, the ground state arising from 
spontaneous symmetry breaking is also characterized by 
a phase Q . 

Traditionally, lasers and superconductors are studied 
separately. Recently 0], it has been realized that the 
superconducting (SC) phase difference and the optical 
phase may interact in a single device that combines two 
superconductors and a semiconducting p-n junction. The 
latter is a common system for light generation as the 
electron-hole recombination produces photons of visible 
frequency Combining semi- and superconductors 

within a nanostructure has been a difficult technologi- 
cal problem that attracted attention for a long time [5|. 
It has been solved using semiconductor nanowires @ or 
quantum wells Q, opening up the possibility to make 
combined devices. 

The device in question has been termed a Josephson 
LED [Fig.QJa)]. It employs a double quantum dot (QD) 
in a p-n semiconductor nanowire connected to SC leads 
Q. The device, biased with a voltage V, exhibits two 
types of photon emission: "blue" photons at the Joseph- 
son frequency wj = 2eV/h due to the recombination of 
a Cooper pair from each side of the junction, and "red" 
photons at about ujj/2 due to electron-hole recombina- 
tion. It has been shown that the optical phase of the 
Josephson generated "blue" photons is locked with the 
SC phase difference. The resulting "blue" light could in 
principle be enhanced by traditional optical methods but 
its small intensity makes this a challenging task. 

In this Letter, we explore an alternative idea where the 
far more intense "red" emission is enhanced in a resonant 
cavity mode. We find lasing at half the Josephson fre- 
quency and, thus, dub the device 'Half- Josephson Laser' 
(HJL). In a common laser, lasing results from sponta- 
neous symmetry breaking where all values of the optical 
phase are equivalent. Drift between these values leads to 



a finite decoherence time. In contrast, the optical phase 
of the HJL is locked to the SC phase difference with only 
two allowed values of the optical phase corresponding to 
two opposite radiation amplitudes. This removes drift as 
a source of decoherence and opens up the possibility to 
manipulate the optical phase by changing the SC phase 
difference. Instead, decoherence of the radiation in the 
HJL results from switching between different QD states 
accompanied by the emission of a photon. We have ex- 
plored these processes and find that by order of magni- 
tude the resulting decoherence time is the same as the 
theoretical limit for a common laser Td oc = n/T, with T 
the damping rate and n the number of photons accumu- 
lated in the resonant mode. A rather low T is required to 
achieve lasing for a single Josephson LED, this condition 
being relaxed with a large number of LEDs in a single 
cavity [§]. 

Setup and model The HJL is a Josephson LED em- 
bedded in a single mode optical cavity with resonance 




FIG. 1. (Color online) (a) The Josephson LED: electron 
(above) and hole (below) QD levels are close to the chemical 
potentials fi Ey h of the SC leads which differ by an energy eV. 
Charge transfer is only possible either through electron-hole 
recombination with the emission of a "red" photon at icjj or 
through a Cooper pair transfer with the emission of a "blue" 
photon at uij. (b) The HJL is a Josephson LED embedded in 
an optical cavity with a resonance frequency uuo « eV/h, i.e., 
close to the "red" emission frequency. The separately-colored 
regions in between the depleted areas represent the two QDs. 
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frequency u>o ss Fig. [TJb). The light emission from 
the cavity is described by a damping rate T. The elec- 
tronic part consists of a biased p-n junction where each 
side of the junction accommodates a QD connected to 
a SC lead. The barriers separating the QDs from the 
leads are arranged such as to allow charge transfer only 
through electron-hole recombination. Such QD junctions 
can be realized with semiconducting nanowires [9j. 

The minimal model for the QDs involves a single or- 
bital for each QD. An orbital can house up to two parti- 
cles (including spin) yielding 16 possible states. The QD 
Hamiltonian then reads [l(J 

Hqo = ^ \. E &i + Uihi{hi - 1)] + U eh n e n hl (1) 



hilation) operator ¥ (b). Thereby, the photon-dependent 
part of the Hamiltonian reads 



i—e,h 



where h e = c a c a (h h = ^IK) is tne electron 
(hole) number operator, and c a (h a ) is the annihilation 
operator for an electron (hole) with spin a. The energies 
E ej h are measured with respect to chemical potentials 
jU e ,ft of the corresponding leads that differ by an energy 
eV = fx e — Hh', here, U e ,h > is the on-site charging 
energy and U e h < the Coulomb attraction between 
electrons and holes. For concreteness, we assume that the 
hole level houses a h eavy hole with J z = ±|?i, where z is 
the nanowire axis [J, . Such levels are commonly used 



in optical experiments with QDs 12|. Our qualitative 
results do not depend on this particular choice. 

Due to the proximity of the SC leads, Cooper pairs 
can coherently tunnel between the SC leads and the 
QDs introducing mixing between unoccupied and dou- 
bly occupied QD states. These processes can be com- 
pactly described by an additional term H sc = A*c^c± + 
Ahhfhi + H.c. in the Hamiltonian; here, the induced 
pair potentials A e< h have reduced magnitudes in compar- 
ison with the gaps A e) /j of the SC leads but they retain 
the same phases <fi e ,h- Owing to gauge invariance, the 
physical quantities depend only on SC phase difference 
<j> = 4> e — (f>h The Hamiltonian is valid under the condi- 
tions |A e ,ft|,-E e ,fc! U e ,h, U e h ^5 We note further that 
this Hamiltonian along with the electron- hole recombina- 
tion conserves parity (even or odd) of the total number 
of particles on the QDs. Even-odd transitions require 
creation of quasiparticle excitations in the SC leads and 
occur with a relatively slow rate estimated below. 

Interaction between the resonant mode and QDs is de- 
scribed by -ffi n t = — E • d; E being electric field of the 
mode at QD position and d the dipole moment of the 
optical transition between the conduction and the va- 
lence band. We assume a linear polarized mode, choose 
the x-axis in the direction of the polarization, and no- 
tice that for heavy holes d x oc (xe~ %eVt l h + H.c.) with 
x = (h^df + hfti). The time-dependence of the dipole 
moment is due to the applied voltage. It is convenient 
to implement a rotating- wave approximation transferring 
the time-dependent factor to the photon creation (anni- 



(2) 



with ui being the frequency detuning, uj = ujq — eV/h, 
\cj\ <C We see that x plays the role of a driving 

force that excites the oscillations in the mode. We note 
that all Hamiltonians considered conserve spin. 

Semiclassics The present model is a rather complex 
case of nonequilibrium dissipative quantum mechanics. 
However, since we envisage a large number of photons 
in the mode, we employ a semiclassical approximation 
replacing b i-> (6) = X/G. The Hamiltonian H QD + H SC + 
i/ph can then be diagonalized to obtain the spectrum 
E m (X) and corresponding eigenstates \m). The dipole 
strength x m ( A) = (m\x\m) depends both on the radiation 
field A and the QD state \m). Since the dipole strength 
in turn determines the evolution of the radiation field via 
the evolution equation 



A = - 



A 



G 2 

i~ x m (A) , 



dE„ 



(3) 



we have to solve the system self-consistently [l[ . The ra- 
diation field can build up as long as the energy gain rate 
2HuJo(G 2 /h)Im[x m (X)/X\ due to the nanowire is greater 
than the energy loss rate ftwoT. With increasing A the 
energy gain saturates till a stationary state of radiation 
(SSR) with A = is reached at a certain radiation am- 
plitude A s . 

In conventional lasers, the driving is due to a popula- 
tion inversion that originates from dissipative transitions 
in an open system. For the HJL, the SC drive is not 
dissipative by itself: only the emission of photons from 
the cavity is a dissipative process. The driving origi- 
nates from coherent mixing of discrete quantum states 
due to the proximity of the QD to the SC leads without 
any population inversion. Thus, the driving mechanism 
of the HJL is very different from that of a conventional 
laser. This is why the information about the SC phase 
difference is preserved in the process of driving. The en- 
ergy gain, including its sign, depends on the difference 
between <f> and the phase of A. Owing to this, the phase 
of A s of the SSR is locked to the SC phase difference. The 
SSRs of the HJL come in pairs ±A S which is very differ- 
ent from a conventional laser where only the magnitude 
| A s | (photon number) is fixed. We give in [§] analytical 
solutions to Eq. ([3]) for a toy two-level model. 

Scales Let us estimate the scales involved that are ex- 
pected to yield lasing. To simplify, we assume all char- 
acteristics of the QD spectrum to be of the same en- 
ergy scale E which is of the order E ~ \&e,h\ <^ 
This assures optimal mixing of the QD states by super- 
conductivity. In a lasing state, the radiation amplitude 
should noticeably contribute to the energies of the QD 
states. This requires |A| ~ E. Assuming uj ~ T we esti- 
mate from Eq. ([3]) that this takes place at G ~ \fhTE. 
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We will assume that G is always chosen to be of this 
scale. The number of photons is then estimated as 
n ~ \X\ 2 /G 2 ~ E/hF. The semiclassical approximation 
is thus justified provided T is sufficiently small, T -C E/H. 

Lasing Despite the model being minimal, it contains 
ten parameters that affect the existence and characteris- 
tics of the SSRs. To find these characteristics, we need 
to evaluate the dipole moment in a given state at given 
A, and can do it separately for the states of odd and even 
parity since they are not mixed by interactions. Addi- 
tionally, the spin conservation splits the eight odd states 
into two equivalent groups of four corresponding to total 
spin ±i. For the even states, only one of the four pos- 
sible |l e l?i) states, + Mct)|0), couples to the field. 
Hence, we only need to consider five of the eight even 
states as the other three are dark. 

With this, we demonstrate lasing as proof of concept by 
finding SSRs in the even states for QD parameters within 
the above estimated scales, — E e = E) x = \U e — \Uh = 
—U e h = A/j = E and A e = 1.5E, for wide regions in the 
space of detuning oj and coupling G, see Fig. [2j Note 
that each eigenstate |m) has a different dipole strength 
x m such that the lasing threshold G c [Fig. EJa)] and the 
radiation amplitudes of the SSRs A s [Fig. [^Jc) and (d)] 
depend on m. Figure [^fb) shows the number of photons 
upon crossing the lasing threshold for the state |2). In 
agreement with the estimations, n reaches the maximum 
~ E/HT at G ~ vhFE. We stress that the optical phase 
of the radiation amplitude in an SSR is not arbitrary but 
locked to the SC phase difference [Fig. (He)]. 

Switching In the above discussion, we have assumed 
the QD to stay in a certain eigenstate \m). In fact it 
does not: the finite spectral width of the mode enables 
switching between the eigenstates. As shown below, the 
switching events occur on a much longer timescale 
than that of the relaxation of A towards its stationary 
value, r _1 . This separation of timescales allows to con- 
sider the switching dynamics separately from the dynam- 
ics of the radiation amplitude. 

Each switching event is accompanied by the emission 
of a photon with a frequency mismatch compensating the 
difference of energies between initial and final eigenstates, 
fiWk = Ef — Ei. For switchings not altering the parity 
of the eigenstate, the rates r sw can be evaluated using 
Fermi's Golden Rule that contains the effective density 
of photon states T/Huj 2 (the tail of a Lorentzian-shaped 
emission line of the resonant mode) and the square of 
the matrix element, |(m/|i? p h|rrai)| 2 , with \mnj)) denot- 
ing the initial (final) state. Thereby, the switching rate 
can be estimated as r sw ~ TG 2 /E 2 ~ T/n -C T. The 
switching events are thus rare and the device stays in one 
of the SSRs between the events. 

Switchings altering parity are even rarer as they re- 
quire the excitation of a quasi-particle above the SC en- 
ergy gap lAe^l- The larger detuning of the off-resonant 
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FIG. 2. (Color online) SSRs in the even states for QD param- 
eters given in the text. The five eigenstates are labeled with 
numbers, (a) Lasing thresholds G c for three eigenstates. At 
the line with the number m the energy gain at A = for the 
state \m) exactly equals the energy loss, (b) Number of pho- 
tons n for eigenstate |2) versus the coupling constant G for 
ui/Y = —0.1 [dotted line in (a)] above the lasing threshold at 
G = G c ~ 0.5v hTE. Plots (c) and (d) illustrate the radiation 
amplitudes A™ (marked with circles) for SSRs corresponding 
to the different eigenstates of the QD. The parameter choice 
is given by the cross in (a) where only the states |2) and |3) 
are lasing. (c) Non-zero A™ come in pairs with opposite sign. 
Changing the SC phase difference will rotate the A™ with re- 
spect to the origin of the plot, (d) Eigenenergies versus A (at 
ImA = 0). The A™ are different for each \m). 



photon uik — \&-e,h\/h and an additional small fac- 
tor |A e ^/A ej ^| result in a paramctrically smaller rate 
r _ ~ |A|rG 2 /|A| 3 < T sw [3]. Such processes do not 
conserve spin thus enabling switchings between dark and 
emitting states. 

It is important to realize that, since the SSRs for dif- 
ferent eigenstates have different values of A™, A docs 
not jump to the new stationary value upon a switch- 
ing. Rather, the amplitude will evolve to A s f within a 
timescale ~ , according to Eq. ([3]). For the same 
reason a switching event always involves different eigen- 
states rather than different SSRs at the same eigenstate. 
The latter would require large fluctuations of A that are 
suppressed exponentially. Fig. [3] shows a sketch of the 
radiation intensity as function of time. In contrast to 
common lasers, the HJL intensity fluctuations are large 
at timescales of T~^. 

Decoherence The intrinsic mechanism of decoherence 
in common lasers is a drift of the optical phase. For the 
HJL, this mechanism does not work since the amplitudes 
of the SSRs are locked to the SC phase difference. This 
renders switching the most important source of decoher- 
ence in the HJL. Indeed, after switching from a lasing 
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FIG. 3. (Color online) Sketch of the radiation intensity |A| 2 
evolving in time, (a) Switching events not altering the parity 
of the QDs occur at the timescale ~ F^l. After switching, 
the radiation amplitude attains its new stationary value at 
a timescale ~ <C ^svv during which the QD remains in 
the same eigenstate. (b) Switching events altering the parity 
occur at a longer timescale ~ F"* S> ^sw- These can change 
between the dark and lasing states. 



to a non-lasing SSR the radiation extinguishes quickly 
and its phase is forgotten. Even if the next switching 
brings the system to a lasing eigenstate, the radiation 
will evolve from the initial A = to any of the two pos- 
sible ±A™, with equal probability. Since decoherence 
is due to switching, the relevant timescale is given by 
Tdcc — ~ n/T. Despite the very different decoher- 



ence mechanism 
common laser 13 1. 



this estimation is the same as for the 



Average power and current The intensity fluctuations 
due to switching self-average at timescale exceeding T~} o . 
The averaged characteristics are expressed in terms of 
the probabilities PJ™ to be in a SSR s that belongs to an 
eigenstate |m). Those are given by the stationary solu- 
tion to the master equation of the switching dynamics, 
that is composed of the switching rates Q. In terms of 
these probabilities, the average number of photons in the 
cavity is given by n = J2 m s ET'l^Tl 2 /*-* 2 - The average 
emission power is proportional to the photon number, 
W = TjhuijTn. The same holds for current in the de- 
vice: since the emission of each photon is accompanied 
by a charge transfer, it is given by I = eTh — W/V. 
An elaborated example of the current/intensity depen- 
dences is provided in [§]. For the current-voltage char- 
acteristic, we find a rather complex structure beside a 
peak with a magnitude of the order eE/h, that is con- 
centrated in a narrow interval ~ HT/e of voltages in the 
vicinity eV = hujQ. In this structure, two types of discon- 
tinuities are present: (i) kinks marking the thresholds of 
lasing instability at A = (second-order transitions) and 
(ii) jumps signaling the appearance of a SSR with sta- 
tionary radiation amplitude A s far from A = (first-order 
transitions). We observe a relatively high probability to 
remain in an SSR with a large photon number. It is 



explained from the fact that eigenstates at large A are 
close to eigenstates of x. This suppresses off-diagonal 
dipole-matrix elements resulting in a suppressed rate of 
transitions from this state. 

Feasibility To show the feasibility of the HJL, we 
present here the estimations with concrete numbers. The 
SC gaps | A e> h| are typically ~ 1 meV, so we can choose 
the QD energy scale E ~ 0.1 meV. To estimate the dipole 
strength G ~ ea|E | , with |E | ~ \J hu a /Vol being the 
quantum fluctuation of the electric field in the mode, we 
assume the cavity volume Vol ~ £ 3 with the wavelength 
I = 2ttc/ojq ~ 600 nm, and take a ~ 5 A for the atomic 
distance scale. This gives the maximum G ~ 0.1 meV. 

With these two values for E and G, the minimum 
damping rate required for lasing is T ~ G 2 /hE ~ 
10 11 Hz, corresponding to quality factor Q ~ 10 3 which 
is common for optical cavities. However, in this situ- 
ation the number of photons n ~ 1. This can be en- 
hanced by increasing Q and simultaneous decreasin g G 
so it remains ~ s/HTE. For photonic crystal cavities [14[ 
quality factors Q ~ 10 6 [iff have been measured and 
Q ~ 10 s (HI have been theoretically predicted. This 
gives photon numbers n ~ 10 3 and n ~ 10 5 , respectively. 
The estimations of the emitted power and current at the 
peak do not depend on the choice of T and are given by 
W ~ 10 nW, I ~ 10 nA. The requirements on T can be 
eased and n enhanced by putting many Josephson LEDs 
in the same cavity. Furthermore, this also increases the 
emission power W and the current I 

In conclusion, we have demonstrated the feasibility of 
generating coherent visible light at half the Josephson fre- 
quency in a SC nanodevice. The workings of the device 
resemble the spontaneous parametric down-conversion in 
nonlinear optics [l7| with the superconductors playing 
the role of coherent optical input. The novel driving 
mechanism results in locking between optical phase and 
SC phase difference. The decoherence of the emitted light 
originates from the switchings between different quantum 
states of the device. 

We acknowledge fruitful discussions with N. Akopian, 
L.P. Kouwenhoven, M. Reimer, and T. van der Sar and 
financial support from the Dutch Science Foundation 
NWO/FOM. 
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Supplementary Material to 'Proposal for an optical laser producing light at half the Josephson frequency' 

This material consists of three parts. In the first part we show explicitly how to find stationary values A™ of the 
SSRs from Eq. <j3j) using a simple toy two-level model of the QD. In the second part, we provide a detailed example of 
the averaged quantities, lasing intensity and current, in dependence on voltage as promised in the main text. Finally, 
we briefly consider a device made by placing many Josephson LEDs into the same optical cavity. 



Toy two-state model 

A presentational problem with the setup described in the main text, is a relatively large number of QD states even 
when using all conservation laws (4 states for the odd parity or 5 states for the even parity). To circumvent this 
problem and still discuss the essential properties of the self-consistency equation and the novel driving mechanism, 
let us consider a toy two-state model. The great advantage of this approach is that all calculations can be performed 
explicitly. 

Two-level system For the toy model, we take one of the states to be |l e l/i), which is unaffected by the supercon- 
ductivity in the leads. The other state is chosen as a superposition of two states with an even number of particles in 
each dot, cos9\0 e 0h) + sin0e l ^|2 e 2/j). These number states are mixed by the SC leads with an angle 8. The phase <j> 
is the SC phase difference. The Hamiltonian H — Hqu + X*x + Xx> in this subspace is of the form 

- ( E cos0A + sin0e^A*\ 

ycos0A*+sin0e-^A -E )' ^> 

Here, 2E is the energy difference between the two states in the absence of radiation and A = G(b) is due to the dipole 
coupling to the resonant mode. Diagonalizing this Hamiltonian yields the eigenenergies 



E±(X) = ±\jE 2 + \X\ 2 + sin(20)Re[A 2 e-^], (5) 

where '±' labels the two eigenstates of Eq. (QJ. Note that for A = we have |+) = |l e l/i) and |— ) = cos0|O e Oft) + 
sin6'e l *|2 e 2ft). The average values of the dipole operator in the eigenstates are calculated as 

*±(A) = <±|*|±> = ^ = ± , A + sil ^ A * (6) 
dX* 2^E 2 + [Ap + sin(2(9)Rc[A 2 e-^] 

The dipole serves as a driving force in the self-consistency equation which for stationary states is given by 

r\ c 2 

,+ -)x = - l — X± (X). (7) 

As was noted in the main text, a laser field can develop when the energy gain rate Ge = 2hwo(G 2 / h)lm[x m (X) / X] 
due to the Josephson LED is larger than the energy loss rate HujqT. With increasing A the energy gain will saturate 
because the denominator of Eq. ^ increases, until a stationary state of radiation (SSR), A s , is reached (A s = 0) and 
steady state lasing is achieved. For this toy model the energy gain is given by 

_ G 2 sin(2g)sin(2 X -0) 
E± h ^E 2 + \X\ 2 [1 + sin(20) cos(2 X - 0)] 

where x is the optical phase, A = |A|e Ix . 

Driving mechanism In the main text we have stressed the novelty of the HJL driving mechanism. In conventional 
lasers, the driving is due to a population inversion that originates from dissipative transitions in an open system. 
However, for the HJL the drive not dissipative and the energy gain depends crucially on mixing of QD states by the 
induced SC gaps without any population inversion. 

In our toy-model the essence of the driving comes about through the mixed state cos0|O e O/j) + sin 9 e l ^\2 e 2h) which 
renders both A*(— |x|+) and A*(+|x|— ) non-zero. Therefore, irrespective of the eigenstate of the QD, the system can 
always decay to the other state and emit a photon. Without mixing (0 = 0, it) this is not possible. The energy required 
for the decay is supplied by the bias. With increasing A the values of the amplitudes decrease as the eigenstates | + (— )) 
will become increasingly more like the eigenstates of x. Hence the driving saturates. 
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From the expression for the energy gain, we explicitly see the role of the SC phase difference, <j>. First, the value of 
this phase determines the sign of the energy gain, and thus whether it acts as a gain or as an absorber of radiation. 
Second, the SSRs with radiation amplitudes, A s , depend on the specific combination of phases, 2\ — <f>- Hence, the 
SSR values of the optical phase must depend on the SC phase difference, which signifies a phase-lock between the 
optical phase and SC phase difference. 

We note that, because of this specific combination of phases, there are always two optical phases possible since 
changing x — > x + 7r yields an SSR with the same radiation magnitude |A S | but opposite sign. Hence, the SSRs come 
in pairs ±A S . This is in stark contrast to conventional lasers where the driving is independent of the optical phase 
and all phases occur with equal probability. 

As a final remark concerning the drive we note that the dipole x± completely saturates in the limit of A 3> E, to 

f + sin(260e^~ 2lx 

x± w ±e x 



2y/l + sin(20) cos(2x - <j)) 

In this limit the eigenstates of H are the eigenstates of x so that the amplitudes (— \x\-\-) and (+|x|— ) vanish. 

Stationary states of radiation To find the radiation amplitudes of the SSRs, A™ =± , we need to solve the stationary 
self-consistency equation with the dipole of Eq. ©. Here we see by inspection that A = is always a stationary 
solution since the dipole is proportional to |A|. It may be unstable however against small fluctuations in the mode. 
By performing stability analysis we distinguish between the two cases. 

To find a SSR with nonvanishing radiation amplitude, we first divide both sides of Eq. Q (taking dipole x+) with 
its conjugate to derive expressions for the optical phase 

2l7 = I + 2ioj/T = I + sin(2g)e^- 2t x 
6 l-2ius/T 1 + sin(20)e 2i x-^ ' [ ' 

where 7 = arctan(2w/r). This reduces to 

COS T 

cos(2x-0 + 7) = --7^. (10) 

with extra condition sgn[sin(2f?) sin(2x — 4>)] — 1. From this, we explicitly see that the optical phase \ may assume 
two equivalent values shifted by it and that it is locked to SC phase <j). Nontrivial solutions for the optical phase x 
exist for us < — ir|cot 20 1, i.e., only at negative detuning. For the '— ' eigenvalue, the situation is reversed in detuning: 
the lasing solutions exist only at positive us with us > ^T\cot29\. 

Now taking modulus square of both sides of (JTJ) , we arrive at an expression for the optical field magnitude 

G 4 _ 4 ff 2 + |A| 2 [I + sin(2fl)cos(2 X -0)] 



h 2 (uj 2 + |r 2 ) 1 + sin 2 (26>) + 2 sin(20) cos(2x - 0)) 

Solving Eqs. (jTOj) and (JTTJl simultaneously yields all SSRs, A™ =± . 

From Eq. (jlll) we may derive expressions for the lasing thresholds as it implies that the lasing solutions with | A | > 
can only exist above some critical coupling 



G| 2|£|^ 2 + ir 2 m 

n 



yjl + sin 2 (2(9) + 2 sin(20) cos(2x - 4>) co&2 ^ 



u>\ - Jsm 2 (26)uj 2 - \ cos 2 (26»)r 2 



Stability analysis shows that the zero stationary solution becomes unstable at this critical coupling. Because the lasing 
field continuously starts to grow when crossing this threshold, we have the analog of a second-order phase transition. 
Interestingly, at yet higher coupling 



G 2 c2 2\E\ 



cos 2 (26») 



sin 2 (26V 2 " |cos 2 (26>)r 2 



the zero solution becomes stable again and we encounter an analog of the first-order transition where the SSRs with 
nonvanishing radiation amplitudes still exist at the threshold. At the tricritical point 



the transitions of both orders coexist. 



8 



even parity 



G = JhYE 



odd parity 




0.04 



53 



/ 4 




(b) ; 


\ 










/ 


\ 







-i co/r l 



-2 -1 co/r 1 2 



FIG. 4. (Color Online) Time-averaged current versus detuning for the even and the odd parity states. The dashed curves 
(reduced by factor of 3) give the current in individual lasing QD eigenstates, with the numbers in the even parity plot matching 
those of Fig. [5] in the main text. The solid curves represent the time-averaged currents which are the individual currents 
weighted with the probabilities P s m , the stationary solutions to master equation (|12[) of the switching dynamics as explained 



in the text. The QD parameters are as given in the main text, 
and G = VhTE. 



-E e = Eh 
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A h = E and A e = 1.5E, 



Example of current /intensity dependence on voltage 

In the main text the dynamics of the HJL over longer times was described in terms of switching rates. It was argued 
that these switchings occur after such long times that the corresponding dynamics is decoupled from the dynamics of 
the lasing field A. Hence, the switchings always occur when the system is at an SSR of Eq. ([3]), A™, where the QD 
is in eigenstate \m). We additionally argued that (for |A| -C |A|) switchings that change the parity of the QD state 
occur at an even slower rate. Even for |A| ~ 0.1|A| the parity switchings r c _ D are an order of 10 3 slower than the 
switchings T sw . Therefore, when measuring over times much longer than T~-^ but shorter than T~^, we see an average 
intensity corresponding to the SSRs of many QD eigenstates with equal parity. After a typical time scale T~l the 
QD states switch parity and the average intensity corresponding to this parity is seen. The purpose of this section is 
to calculate the average intensity as a function of voltage for the even and odd parities separately and discuss some 
important features found in the corresponding curves. 

The average intensity can be found by solving the master equation governing the switching dynamics and finding 
the corresponding stationary solution. The probability p™ to be at A" 1 evolves in time according to 

X) W »».WP? (12) 

where W m ,s;n,r denotes the transition from A™ t— > A™ and we define the transition matrix W. The summation is over 
the stationary solutions A r of all eigenstates, \n). To solve master equation (|12|) we need to find and diagonalize W. 
We notice that this matrix has a left eigenvector with zero eigenvalue because J2 m s ^m,s;n,r = 0. Hence there must 
also be a right eigenvector with zero eigenvalue, which is the stationary distribution of probabilities (p™ = for all 
m) when normalized to unity. Since we are only interested in these stationary solutions we only need to find the right 
null vector of W. We represent this stationary solution as P s m , with a capital P. Hence, we can find P s m when all 
Wm,s;n,r are known. 

To calculate the transitions W m ^- n ,r we note that the switchings are accompanied by the emission of an off-resonant 
photon so that the rates follow from Fermi's Golden Rule. The transition rates are then proportional to the product of 
\(m(\?)\H ph \n{\?))\ 2 and a Lorentzian shaped photon density of states, p(AE) = (l/2ir)hr/[(hT /2) 2 + (AE) 2 }, with 
AE = \E n — E m \ the energy difference between |n(A™)) and |m(A")). We note that the transitions go from eigenstate 
|n(A™)) H> |m(A")) (n ^ m), which are the eigenstates of the Hamiltonian _ff QD + H sc + H p h in the semiclassical 
approximation, b t— > (b) = X/G. Additionally we stress that the value of A is unaltered during the transition, which is 
indicated by the (A™) in (m(A")| and |n(A™)). Only after the switch will A evolve to A™ according to Eq. (|3]). 

To find the matrix elements of H p h [Eq. ([2])] we only need to consider the interactions G(Wx + bx^) as the term 
hujwb will always yield zero contribution when n/m. The two interactions imply two possible transition amplitudes 
distinguishing two types of emission: one far below and one far above the cavity resonance. Since we are in the 



dt 



E 

n,r^rn,s 



[W„ 
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rotating frame, this translates to creating a negative frequency photon with bx\ or creating a positive frequency 
photon with wx. We do not expect absorption of off-resonant photons as they are very scarce. Taking the detuning 
into account, the photon energies are tkJ^ = fku — AE and fkuk = Hu) + AE respectively. The detuning may be 
neglected, however, as it is typically taken of the same order as T -C AE. For this same reason we may also neglect 
the r in the denominator of the density of states p(AE). 

Taking the above considerations into account we find the transition rates 

VG 2 

W m ,s;n,r = j^\{m{\ n r )\x\n{X n r )) (m(A™)|^| n ( A ™))| 2 . (13) 

With this we can find W and thus P"\ We make some brief remarks about W. As stipulated above the switchings 
only occur between different QD eigenstates. Fluctuations that let the system jump from one to another SSR in 
the same QD eigenstate are exponentially small due to the large number of photons. Thus W mjS;m r = if s ^ r. 
Furthermore, we note that when in a lasing eigenstate and starting at A = the system may evolve to both SSRs 
with ±A™ with equal probability. Hence the transition rate from a non-lasing QD eigenstate to the SSR with A™ of 
a lasing QD eigenstate is half the value calculated in Eq. (fT3| . For transitions between two lasing eigenstates it must 
be determined to which of the multiple stationary points A™ the field evolves. This is done by integrating the time 
evolution equation of A, Eq. ([3]), until one of the A™ is sufficiently approached, using the radiation amplitude A" of 
the previous SSR as initial condition. 

With the stationary probability distribution, P s m , we calculate some measurable quantities. First, the average 
number of photons in the cavity is n — J2 m s -PJ" | A^ | 2 / G 2 . Then the average number of photons that escape the 
cavity is Tn, yielding the average emission power W = ^HujjTn and the average current in the device / = eTn. The 
latter because every photon emission is accompanied by a single charge transfer. 

Using the minimal model of the main text we calculate the current /intensity as a function of detuning. As mentioned 
above we plot the dependencies for odd and even parity states separately. The results are presented in Fig. 21 Since 
the detuning is related to the applied voltage via uj = eV/h — uq, the peaks in the current, of the order of eE/H, are 
in a narrow interval ~ hT/e of voltages. There are two kinds of discontinuities in the curves: kinks, that mark the 
thresholds of the lasing instability at A = in some QD eigenstate (second-order transitions), and jumps, that signal 
the appearance of an SSR far from A = (first-order transitions). The curves also show a relatively high probability 
to be in lasing QD eigenstates if available. This can be qualitatively understood by considering that eigenstates at 
large A are close to eigenstates of x. This suppresses off-diagonal matrix elements, Eq. (|13p . resulting in a slower rate 
of transitions away from a lasing state. For the odd parity states in Fig. [4J the time-averaged current drops to zero 
while a lasing state is still available. This is because the non-lasing solution in this lasing eigenstate becomes stable 
at |w|/r < 1 and whereas all other states are nonlasing. The device then gets stuck in the SSRs at A = even though 
SSRs with A s ^ are available. 

Many Josephson LEDs in a single cavity 

We have discussed in the main text that the large number of photons n ~ E/hT (assuming G ~ yJhTE) in a cavity 
with a single Josephson LED requires high Q factors. This requirement can be easily softened, and the lasing can be 
achieved at higher damping rate T, if we incorporate N ^> 1 Josephson LEDs in the cavity. Since the LEDs are of ten 
nanometer scale, one can put hundreds of them even into a single-wavelength cavity. The resulting system certainly 
deserves to be explored in detail, and will be a topic of separate research presented elsewhere. In this note we just 
give some straightforward estimations. 

All N QD dipoles are in this setup coupled to a single mode, were we assume equal coupling coefficients G. The 
self-consistency equation for the radiation amplitude changes to 

\= - \-i^Y^x mi (X), x mi = ~Tjjir- (14) 

Assuming for simplicity all dipoles to be in the same direction with same or similar parameter values, we may simplify 
the sum over dipoles to Nx m . This yields A ~ G 2 N/hT for the SSRs. The optimal regime near the lasing threshold 
is achieved if A ~ E. Thus this corresponds to G ~ yj KTE/N. Therefore, the lasing threshold can be achieved at 
much larger damping rate T ~ G 2 N/hE. 

This also increases the number of photons in the mode. If we optimize the coupling to G ~ y/HTE/N, the number 
is estimated asn~ NE/hT, N times larger than the estimate for a single-LED device. If we fix V, G and increase the 
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number of wires in the device, the photon number scales as n ~ N 2 G 2 /T 2 . Note that in both cases the emitted power 
increases as well. It is respectively given by W ~ TeVn ~ eVNE/h, for optimal coupling, and W ~ eVN 2 G 2 /T 
when only increasing the number of wires. 

It is likely that a large number of LEDs in the cavity greatly reduces the intensity fluctuations. With only one 
pair of QDs the lasing intensity fluctuates down to zero as can be seen from Fig. [3J in the main text. In contrast, 
the switching of only one of N QDs will hardly affect the intensity, so the relative intensity fluctuations are expected 
to be of the order of 1/y/N. Likewise, the phase of the optical field will be barely affected by the switching in one 
nanowire. Therefore the decoherence times may become significantly longer than in a single LED HJL. 



